RELATIVISTIC CONIC BEAMS AND SPATIAL 
DISTRIBUTION OF GAMMA-RAY BURSTS 



Heon- Young Chang and Insu Yi 
Korea Institute For Advanced Study, Seoul, 130-012, Korea 

hyc@kias . re .kr , iyi@kias.re.kr 

Received ; accepted 



- 2 - 



ABSTRACT 



We study the statistics of gamma-ray bursts, assuming that gamma-ray bursts 
are cosmological and they are beamed in the form of a conical jet with a large 
bulk Lorentz factor ~ 100. In such a conic beam, the relativistic ejecta may 
have a spatial variation in the bulk Lorentz factor and the density distribution 
of gamma-ray emitting jet material. An apparent luminosity function arises be- 
cause the axis of the cone is randomly oriented with respect to the observer's 
line of sight. The width and the shape of the luminosity function are determined 
by the ratio of the beam opening angle of the conical jet to the inverse of the 
bulk Lorentz factor, when the bulk Lorentz factor and the jet material density 
is uniform on the photon emitting jet surface. We calculate effects of spatial 
variation of the Lorentz factor and the spatial density fluctuations within the 
cone on the luminosity function and the statistics of gamma-ray bursts. In par- 
ticular, we focus on the redshift distribution of the observed gamma-ray bursts. 
The maximum distance to and the average redshift of the gamma-ray bursts are 
strongly affected by the beaming-induced luminosity function. The bursts with 
the angle-dependent Lorentz factor which peaks at the center of the cone have 
substantially higher average gamma-ray burst redshifts. When both the jet ma- 
terial density and the Lorentz factor are inhomogeneous in the conical beam, 
the average redshift of the bursts could be 5 times higher than that of the case 
in which relativistic jet is completely homogeneous and structureless. Even the 
simplest models for the gamma-ray burst jets and their apparent luminosity dis- 
tributions have a significant effect on the redshift distribution of the gamma-ray 
bursts. 



Subject headings: gamma rays:bursts 
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1. INTRODUCTION 

The BATSE experiment on the Compton Gamma Ray Observatory and the study 
of the afterglows (e.g., Piran 1999 and references therein), have established that the 
gamma-ray bursts (GRBs) are cosmological (Mao and Paczyhski 1992; Meegan et al. 
1992; Piran 1992). Even though the distance scale seems settled (Metzger et al. 1997), it 
appears that uncertainties remain in the total energy and the burst rate of GRBs (Kumar 
1999; Kumar and Piran 1999). These two important issues depend on the level of beaming 
in GRB emission. That is, the issues critically depend on whether the geometry of the 
gamma-ray emitting ejecta is spherical or jet-like (Harrison et al. 1999; Kulkarni et al. 1999; 
Meszaros and Rees 1999; Sari et al. 1999). A number of authors have studied energetics 
and geometry of the ejecta (Mao and Yi 1994; Rhoads 1997; Panaitescu and Meszaros 
1998; Rhoads 1999; Moderski et al. 2000). It is also important to put constraints on the 
width of the luminosity function by comparing the observed intensity distribution with 
those predicted by a physical model (Mao and Yi 1994; Yi 1994). In essence, the rate, the 
energy, and the luminosity function of GRBs are all closely related to whether or not the 
geometry of the ejecta is spherical. 

Two most frequently quoted statistics in GRB observations are < V/V max > and 
logiV(> F) — logF, where F refers to the peak flux (or peak count rate) and N denotes 
the number of GRBs with fluxes higher than F (e.g., Yi 1994). These two quantities 
contain information on the lumonosity function of GRBs and the spatial number density 
of the sources. A value of < V/V max > consistent with that of an observed sample is a 
necessary condition but not a sufficient condition for a luminosity function which is 

neither directly observable nor theoretically well undertood. The luminosity function of 
GRBs can be obtained for an assumed source distribution n(z) such that the calculated 
logiV(> F) — \ogF fits the observed distribution, and vice versa. The density n(z) refers 
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to the rate of GRBs per unit time per unit comoving cosmological volume. However, due to 
the very nature of N(> F), which is the convolution of n(z) and one almost always 

obtains n(z) for a given such that the theoretical logiV(> F) — logF curve fits the 

observed intensity distribution. Therefore, in order to extract information concerning n(z) 
or $(X), one has to assume one of these two functions or to develop a techinique to separate 
the effects of these two unknown functions (Horack and Emslie 1994; Horack et al. 1994; 
Ulmer et al. 1995; Ulmer and Wijers 1995). It is therefore of great interest to construct 

on the basis of a physical model, which is one of our major goals in this Letter. Since 
there remain uncertainties in GRB engine models, we focus on the consequences of the 
conical beaming without specifying how a beam is formed in a physical engine model. 

Using the first BATSE catalog of gamma-ray bursts (Fishman et al. 1994), Mao and Yi 
(1994) studied the effects of the relativistic bulk motion in a conical beam on the statistics 
of gamma-ray bursts. They found that the luminosity function is naturally introduced by 
the random distribution of the space orientation of the cone axis and that the case of the 
standard candle is not easily distinguished from that of the beaming-induced luminosity 
function with a sharp peak. This is especially the case for large beam opening angle and 
the large Lorentz factor 7, as one may expect. Different Lorentz factors and opening angles 
however result in non-trivial changes for the distances to GRBs and especially the highest 
redshift of or the maximum distance to the most distant GRB for a given sample. For 
instance, the maximum redshift z max increases as the ratio of the opening angle to 7 _1 
decreases. We modify the conical beam model by allowing a spatial variation of 7 and the 
density profile of gamma-ray emitting electrons on the photon-emitting surface of the cone. 

From numerical simulations of relativistic jets (Marti et al. 1997; Renaud and Henri 
1998; Rosen et al. 1999) and observations of the astrophysical jets (Zensus 1997; Spruit 
2000), it is clear that jets do have some significant structure in them and the bulk Lorentz 



- 5 - 



factors evolve as the jets propagate. Therefore, it is plausible to extend the simplest jet 
model such as that of Mao and Yi (1994). In a more realistic jet model, the bulk Lorentz 
factor has a spatial profile at the surface where the observed gamma-ray emission occurs 
and the spatial electron density distribution is significantly inhomogeneous. 

In § 2 we begin with a brief presention of data we use, which are parts of the BATSE 
4B catalog (Paciesas et al. 1999), and we describe our conical jet geometry, following Mao 
and Yi (1994), in § 3 we present results. Finally, we conclude with summary of our results 
and discussions in § 4. 

2. OBSERVATIONAL DATA AND MODEL 

The BATSE 4B catalog (Paciesas et al. 1999) provides 1637 triggered GRBs detected 
from 1991 April through 1996 August. We use the bursts which are detected on the 1024 
ms trigger time scale. We choose the bursts of which peak count rates are above 0.4 
photons cm^s -1 in order to avoid the threshold effects (cf. Mao and Yi 1994). Of those 
bursts, we select the GRBs whose C max /C min is greater than 1.0, which gives a sample of 
651 bursts. The energy range in which the peak flux is measured is 50 — 300 keV. The peak 
fluxes of bursts vary by about 2 orders of magnitude in this data set. 

Throughout this Letter we adopt the simplest cosmological model; the universe is 
flat with the density parameter f2 = 1, the Hubble constant is 65 kms~ 1 Mpc~ 1 , the 
cosmological constant is absent, all the GRBs are 'standard bursts' with an identical 
power-law spectrum, and the rate of bursts per unit comoving volume per unit comoving 
time is constant (cf. Mao and Paczyhski 1992; Yi 1994). 

In our beaming model, the ejecta is flowing outward relativistically in a cone with the 
geometrical opening angle A9. The observed gamma-ray emission is produced at radius R 
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from the central engine where the radiation first becomes optically thin. In the minimal 
beaming model of Mao and Yi (1994), the ejected material has the same bulk Lorentz factor 
7 at this distance and the photon-emitting electrons' density at the surface of the cone is 
uniform. In this Letter those simplifications are relaxed as described below. 

In the cylindrical symmetry one may consider the colatitude alone, which is defined by 
the angle between the line of sight and the symmetry axis of the cone. Besides the spherical 
coordinate system centered on the central engine, we introduce an auxiliary spherical 
coordinate system with the z'-axis along the symmetry axis of the cone. It can be shown 
that the angle 9 between a position of a direction within the cone and the line of sight is 
given by 



The monochromatic flux received by a local observer at a distance D from the source, 
taking the cosmological redshift effects into account, reads 



where Dl(z) = 2c/Hq[1 + z — (1 + -z) 1 ^ 2 ], c and H$ being the speed of light and the Hubble 
constant respectively, and Y = [7(1 — /^cosO)] -1 , (3 — (1 — 7" 2 ) 1 / 2 . And we have used 
the relation v = IV ; here u is the corresponding frequency in the comoving frame. In 
this model we have ignored the structure of light curves due to the relative time delay of 
radiation from different parts of the cone, and other complicated cosmological effects. 

We define the terms, excluding terms for the cosmological information in the above 
equation (2), as the local peak count rate, 



cos 6 = cos 9' cos 9 — sin 0' sin 4> sin 9. 





(2) 
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x / sin#'cos0r 2+c W. (3) 
Jo 

The maximum local peak count rate is achieved at 9 = and can be obtained analytically : 



where x u 1 = 1 — (3 cos A9, and x l 1 = 1 — (3. 

Based on the randomness of the direction of the cone axis with respect to the line of 
sight, we obtain the probability function which directly reflects the angle between the line 
of sight and the direction of the cone, 9. Assuming that the cone is uniformly distributed 
in space, the probability function for a one-sided cone is given by 



where 9 is between and n. Since the local peak count rate is a function of the orientation of 
the cone, one may translate the probability function of the cone's angle into the probability 
function of the local peak count rate. 

In order to obtain the luminosity function and compare with observations, we calculate 
the local peak count rate in a fixed frequency range. We allow variations of 7 and the 
electron density at the surface of the cone by introducing a window function for 7 and 
electron density, N e . The window function is axisymmetric with respect to the symmetry 
axis of the cone. The window function we adopt is the Gaussian function centered at the 
center of the cone: 






(5) 



W(0') = exp[-A(4) 2 ], 



(6) 



where A is a constant, A9 is a given opening angle, and 9' varies from to A9. 
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3. RESULTS 

The probability function of the local peak count rate is shown in Figure 1 for three 
different ratios of A9 to 7 _1 , in the case where 7 is uniform over the photon-emitting 
surface of the cone. For a fixed 7, as the opening angle increases the peak of the probability 
of the local peak count rate becomes higher and narrower, as one may expect. A model 
with a large value of Lorentz factor (A9 ^> 7 _1 ) essentially provides a luminosity function 
indistinguishable from that of the standard candle model. For a given value of the ratio, 
the whole distribution function moves vertically when the absolute value of the Lorentz 
factor changes. This is because the local peak count rate is decreasing faster with the angle 
for larger 7. In the extreme case where A9 <C 7 _1 the distribution is a power law with the 
index of —1/3, as expected from the analytical result (Yi 1993). 

In Figure 2, we show the probability distribution functions for both the varying bulk 
Lorentz factor and the inhomogeneous electron density. We assume the axisymmetry of the 
Lorentz factor profile around the cone axis, 7 = 7(0'), and hence the Lorentz factor profile 
could mimic a simplified model for the jet-environment drag. At the center of the cone, 7 
has the maximum value and decreases with 9'. As 7 decreases more steeply with 9', the 
log Pl (log Pioc) shows a smaller peak at log(p oc / ' Pi oc ,rnax) = 0, and a higher level of the 'tail' 
of the log (log P; oc ). It is because this type of 7 effectively reduces the 'average' value of 
7 over the cone and the 'effective opening angle' simultaneously. For a very narrow window 
function it essentially reduces to the pencil-beam case. The photon-emitting electrons are 
supposed to be distributed according to the 7 distribution such that the local electron 
number density is inversely proportional to the square of the bulk Lorentz factor, 7. That 
is, the electron number density is a function of the angular position in the cone. In this 
profile, the electron number density increases outwards from the center of the cone, since 7 
is a decreasing function. In effect, the cone is reminiscent of the hollow cone. 
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Once the probability distribution is obtained, we are in a position to calculate 
logiV(> F) — logF with an assumed spatial distribution of GRBs, n(z). In Figure 3, we 
compare the cumulative intensity distribution curve produced by our luminosity function 
with observational data. Since the goal of this study is looking at effects of the conic beam 
with varying distributions of 7 and the electron density over the surface of the cone, we 
assume there is no source evolution, that is, n(z) = Hq. We adopt the power law index a of 
unity instead of 2 (cf. Mao and Yi, 1994). This is a simplification in a sense that observed 
GRBs show various power-law indices. This value, however, represents the averaged power 
law indices of observed GRBs (Band et al. 1993; Mallozzi et al. 1996; Preece et al. 2000). 
All the theoretical cumulative probability distribution with the luminosity functions shown 
in Figure 2 are plotted in Figure 3 along with the observed distribution. The curves shown 
in Figure 3 are best-fit functions determined by the Kolmogorov-Smirnov (K-S) test for 
each parameter set. Parameters in our model are not sensitive enough to constrain the 
luminosity function. 

For a given set of model paramters, z max plays a role of a parameter in the K-S test in 
that the luminosity function studied here optimizes logiV(> F) — logF to fit observational 
data. The obtained z max s in this procedure are shown in the Table 1. These z max s obtained 
with the 4B BATSE catalog are greater than those for the IB BATSE catalog by about 0.1 
on average, which indicates that we are seeing fainter GRBs in the 4B BATSE than in the 
IB BATSE catalog and therefore more distant GRBs. 

It shows that the redshift of the most distant GRBs becomes larger as 7 falls more 
steeply from the cone center. As shown in the Table 1, our toy models with the narrow 
Gaussian profiles could easily reproduce z max values as high as the highest reported 
z max = 3.42 (quoted from Bulik (1999), see references therein). This simple case indicates 
that the effects of the luminosity function have significant implications on the cosmological 



spatial distribution of GRBs. It is interesting to see that such a high redshift would be hard 
to explain in the homogeneous beam model with a constant n(z). The effects of the Gaussian 
window on the values of z max are substantial enough for further comments. The model for 
a narrow window with A = 4 gives z max = 3.7(average redshift < z >= 1.03) while the 
uniform beam model (i.e. without any variation of 7) gives z max = 1.6(< z >= 0.51). The 
average value of redshifts, < z >, is taken over the cumulative redshift distribution function 
shown in Figure 4. The cumulative redshift distribution is defined by 



The cumulative redshift distributions shown in Figure 4 also indicate that the 
luminosity distribution induced by the beaming has a significant implication on the GRBs 
distances in the observed samples. When the beam is sharply peaked at the beam center 
(solid curves in the upper panel in Figure 4) the fraction of the high redshift GRBs (e.g. 
z > 3) could be as high as ~ 10% while the broad beam case (e.g. A = 1/8) essentially rules 
out any high redshift GRBs. The cumulative distribution slowly decreases with the redshift 
when the 7 in the conic beam is decreasing rapidly in a sense that GRBs are spead out in 
a broader region beyond the averaged z for the highly concentrated beam. That is, the 
ratio of z max to < z > is 3.59 and 3.14 when A = 4 and A = 1/8, respectively. The beams' 
electron density structure also has a similar effect on the redshifts of GRBs. The effects 
of the luminosity function have to be explicitly considered when observed high redshift 
values are interpreted (Krumholz et al. 1998). For instance, in the standard candle case, 
a significant probability for high redshift GRBs could directly imply a substantial source 
evolution effect. However, the beaming induced luminosity function could make this simple 
interpretation much uncertain (Blain and Natarajan 2000). 



Af(z') 



N(z < z < z max ) 
N(0<z< z max ) 



(7) 



where 
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4. DISCUSSION 

The theoretical models for GRBs are abundant. Despite remarkable progresses in 
understanding physical mechanisms involved in these models, the GRB prompt emission 
mechanisms and engine models have so far been unable to constrain the extent of beaming 
and the luminosity distribution of GRBs. This in turn has been a major uncertainty in 
interpreting the observed flux data in terms of the cosmological spatial distribution of 
the bursts. In this regard, the present work has shown that the simple beaming models 
and their resulting apparent luminosity functions have significant effects in interpreting 
the observed data. If the GRBs are indeed standard candles with a single well-defined 
luminosity, the spatial distribution of GRBs in connection with the cosmic star formation 
rate could be translated into the cosmological source evolution. However, the luminosity 
distributions we have considered affect the maximum redshift and the average redshift 
significantly. It is therefore important to derive a theoretical luminosity function for a given 
GRB model. 

The jet models we adopted are obviously over-simplified. Despite this major drawback, 
the models capture the essential ingredients of the beamed relativistic jets concerning the 
apparent luminosity function. One of the major uncertainties is that the jets and GRB 
sources differ greatly and GRB luminosities and jets' physical conditions are intrinsically 
different in each source. Given the wide range of burst durations and the diverse burst 
types (Fishman and Meegan 1995, and references therein), such a possibility cannot be 
ruled out. If this is indeed the case, our standard source approach is not applicable. 

We thank C. Kim, H. Kim, and K. Kwak for useful discussions and especially C. Kim 
for her help with Figure 3. IY is supported in part by the KRF grant No. 1998-001-D00365. 
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Fig. I. — Logarithm of the probability distribution p L (\ogPi oc ) as a function of logP/ oc . The 
local peak count rate is scaled to its maximum value. In this example, 7 = 10 2 . ^From 
bottom to top, the opening anlges are 0°.1,1°,3°. The photon spectral index a is unity. 
Note that for A9 = 0°.l case, it is basically a power-law function with the index of —1/4. 

Fig. 2. — Plots similar to Fig. 1. In the upper panel, values of 7 varies over the surface 
of the cone. The window function adopted is the Gaussian, exp[— A(-^) 2 ] such that at the 
center of the cone 7 has an original constant value, 7 . Results of this study do not depend 
specific details of the window function. In the lower panel, the electron number density also 
varies over the surface of the cone. The electron density varies as inversely proportional to 
the square of 7 in this example. The thick solid lines correspond to the constant 70 case with 
a — 1, the thin solid lines to A = 2, the dotted lines A = 1, the short dashed lines A = 1/8. 
The opening angle of the cone is 1° and 7 is 100. 

Fig. 3. — The cumulative probability distribution of peak count rate in terms of logarithm 
of the peak count rate. Open triangles are the observed distribution and the solid lines are 
the best fit curves in the Kolmogorov-Smirnov (K-S) test. 

Fig. 4. — The cumulative redshift distribution of bursts with redshifts higher than z as a 
function of redshift z . The upper and the lower panels result from the varying 7 and the 
varying electron density as well, respectively. The solid lines correspond to A = 4, the long 
dashed lines to A = 2, the dot-dashed lines to A = 1, the dotted lines to A = 1/2, and the 
short dashed lines to A = 1/8. 
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Table 1: The maximum redshifts of GRBs which allow the best-fit of models to observed 
data, and the averaged redshifts. The parameter A determines the sharpness of the window 
function. Parameters adopted in the calculation are a — 1, 70 = 100, and A9 = P.O. 





A 




< z > 


constant 70 


N/A 


1.60 


0.51 




4 


3.70 


1.03 




2 


3.17 


0.91 


varying 7 


1 


2.41 


0.73 




1/2 


1.96 


0.61 




1/8 


1.70 


0.54 




4 


11.23 


2.36 
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3.47 


0.98 


varying 7 and N e 
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2.20 


0.67 




1/2 
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